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Abstract 



We study the meson sector of 2+1 dimensional light-front QCD using a 
Bloch effective Hamiltonian in the first non-trivial order. The resulting two 
dimensional integral equation is converted into a matrix equation and solved 
numerically. We investigate the efficiency of Gaussian quadrature in achieving 
the cancellation of linear and logarithmic light-front infrared divergences. The 
vanishing energy denominator problem which leads to severe infrared diver- 
gences in 2+1 dimensions is investigated in detail. Our study indicates that 
in the context of Fock space based effective Hamiltonian methods to tackle 
^^ ■ gauge theories in 2+1 dimensions, approaches like similarity renormalization 

renormalization method may be mandatory due to uncanceled infrared diver- 
gences caused by the vanishing energy denominator problem. We define and 
study numerically a reduced model which is relativistic, free from infrared 
O j' divergences, and exhibits logarithmic confinement. The manifestation and 

violation of rotational symmetry as a function of the coupling are studied 
quantitatively. 

PACS Numbers: ll.10.Ef, ll.10.Kk, 11. 10. St, 11.15.Tk, 12.38. Lg 

I. INTRODUCTION AND MOTIVATION 

There are various well-known motivations |I[ to study QCD in the light-front Hamiltonian 
formalism. In fact, there have been many attempts recently to study relativistic bound state 
problem in the Hamiltonian formalism in a light-front Fock space basis (For a review see, Ref. 
0). It has been realized that a major impediment to a straightforward diagonalization of the 
Hamiltonian is the rapid growth of the dimension of the Hamiltonian matrix with particle 
number. An alternative approach will be to use an effective Hamiltonian that operates in 
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a few particle basis. A challenging problem here is that, for a successful description of 
low energy observables, the effective Hamiltonian must incorporate main features of strong 
interaction dynamics. 

One of the first attempts invoked Tamm-Dancoff || or the Bloch-Horowitz effective 
Hamiltonian. Though it was successful in tackling 1+1 dimensional gauge theories, its defi- 
ciencies become apparent when attempts were made in 3+1 dimensions. First and foremost 
is the lack of confinement in the case of QCD. Second is the appearance of the bound state 
eigenvalue in the energy denominators. This has two undesirable consequences. Firstly, a 
light-front singularity of the type pp, where k + is the light front longitudinal momentum 
of the exchanged gluon, remains in the bound state equation, which would have canceled if 
free energies appeared in the energy denominators. Secondly, from the fermion self energy 
contribution, in addition to the mass divergence another ultraviolet divergence appears (for 
an example in the context of 3+1 dimensional Yukawa model see Ref. Q) which contributes 
to the renormalization of the coupling. This contribution is also infrared divergent and can 
be identified as arising from fermion wave function renormalization. It is the Fock space 
truncation that has produced this unphysical divergence which would otherwise have been 
canceled by vertex renormalization in a strict order by order perturbative calculation. 

It is well-known that various standard formulae for effective Hamiltonian all have draw- 
backs. Some of the deficiencies of the Bloch-Horowitz formalism are absent in the Bloch 
effective Hamiltonian || which was reinvented in the context of renormalization group by 
Wilson |J. Bloch Hamiltonian has two desired properties, namely, the effective Hamilto- 
nian is (1) Hermitian and (2) involves only unperturbed energies in the energy denominator. 
Use of Bloch effective Hamiltonian eliminates two major problems of the Tamm-Dancoff 
approach to gauge theories mentioned above. However, Bloch effective Hamiltonian has an 
undesirable feature, namely, the vanishing energy denominator. To the best of our knowl- 
edge, Bloch effective Hamiltonian was never assessed in terms of its strengths and weaknesses 
in the study of bound state problems in field theory. 

In the study of bound states, QCD poses challenging problems. To overcome many 
pitfalls of standard effective Hamiltonians, similarity renormalization was proposed 0. It 
avoids vanishing energy denominators and thus provides an improvement over Bloch effective 
Hamiltonian. Initial attempts in the similarity renormalization approach worked in either 
the non-relativistic limit || or the heavy quark effective theory context ||. Only recently, 
work has begun [|K| to address many practical problems, especially the numerical ones one 
faces in this approach. 

A major feature of gauge theories on the light-front is severe light-front infrared diver- 
gence of the type , fc K 2 where k + is the exchanged gluon longitudinal momentum which 
appears in instantaneous four-fermion, two-fermion two-gluon, and four-gluon interactions. 
In old-fashioned perturbation theory these divergences are canceled by transverse gluon in- 
teractions. In similarity perturbation theory the cancellation is only partial and singular 
interactions survive. Before embarking on a detailed study of effective Hamiltonian in the 
similarity renormalization approach which is a modification of the Bloch effective Hamilto- 
nian, it is quite instructive to study the Bloch effective Hamiltonian itself. The result of 
such a study can serve as benchmark against which one can evaluate the merits of similarity 
renormalization scheme. This will also provide us quantitative measures on the strengths 
and weaknesses of numerical procedures in handling singular interactions (in the context of 



light-front field theory) on the computer. It is crucial to have such quantitative measures in 
order to study the effects of similarity cutoff factors on the nature of the spectrum. This is 
one of the motivations for the present work. 

Just as the Tamm-Dancoff or the Bloch- Horowitz formalism, Bloch effective Hamiltonian 
of QCD in the first non-trivial order also does not exhibit confinement in 3+1 dimensions. 
Since one of our major concerns is the study of spectra for confining interactions, we go to 
2+1 dimensions. In this case, in the limit of heavy fermion mass, a logarithmic confining 
potential emerges. There are several other reasons also to study light-front QCD in 2+1 
dimensions. They arise from both theoretical and computational issues which we discuss 
next. 

First of all, issues related to ultraviolet divergence become more complicated in the light 
front approach since power counting is different |1[| on the light front. We get products 
of ultraviolet and infrared divergent factors which complicate the renormalization problem. 
Going to two space one time dimensions greatly simplifies this issue due to the absence of 
ultraviolet divergences except in mass corrections. An extra complication is that Fock space 
truncation introduces extra ultraviolet divergences which complicate the situation in non- 
perturbative bound state computations jf|. Such special divergences do not occur in 2+1 
dimensions. A third complication one faces in 3+1 dimensions is that on enlarging the Fock 
space in a bound state calculation, one soon faces the running of the coupling constant. 
At low energy scales, the effective coupling grows resulting in a strongly coupled theory 



TT| making the weak coupling approach with a perturbatively determined Hamiltonian 
unsuitable or making it mandatory to invent mechanisms like non-zero gluon mass to stop 
the drastic growth ]I[. In 2+1 dimensional QCD we do not face this problem since the 
coupling constant is dimensionful in this superrenormalizable field theory and does not run 
due to ultraviolet divergence. We can keep the coupling arbitrarily small and study the 
structure of the bound states in a weakly coupled theory. 

Secondly, in 1+1 dimensions, in the gauge A + = 0, dynamical gluons are absent and 
their effect is felt only through instantaneous interactions between fermions. Further, recall 
that in light front theory, vacuum is trivial. As a result, the Fock space structure of the 
bound states are remarkably simple. For example, the ground state meson is just a qq pair 
both at weak and strong couplings. In contrast, in 2+1 dimensions, one component of the 
gauge field remains dynamical and one can systematically study the effects of dynamical 
gluons. Also note that 2+1 dimensions are the lowest dimensions where glueball states 
are possible and offers an opportunity to study their structure in the Fock space language 
without additional complications of 3+1 dimensions. 

A third reason deals with aspects of rotational symmetry. 2+1 dimensions offer the first 
opportunity to investigate violations of Lorentz invariance introduced by various cutoffs 
(momenta and/or particle number) in the context of bound state calculations. This is to 
be contrasted with 1+1 dimensions where the sole Lorentz generator, namely boost, is kine- 
matical in light-front field theory. Since in 2+1 dimensions we have a superrenormalizable 
field theory, violations introduced by transverse momentum cutoffs are minimal. Thus in 
contrast to 3+1 dimensions, one can study the violations caused by truncation of particle 
number alone and longitudinal momentum cutoffs. It is also conceivable that one can en- 
large the Fock space sector and investigate their effect on restoring Lorentz invariance. It 
is expected that such investigations are more viable in 2+1 dimensions compared to 3+1 



dimensions due to less severe demand on computational resources. 

A fourth reason concerns similarity renormalization approach. In 3+1 dimensions it has 
been shown that similarity renormalization group approach M to effective Hamiltonian in 
QCD leads to logarithmic confining interaction [12|]. It is of interest to investigate corre- 



sponding effective Hamiltonian in 2+1 dimensions especially since the canonical Hamiltonian 
already leads to logarithmic confinement in the nonrelativistic limit in this case. It is also 
known that in 3+1 dimensions the confining part of the effective Hamiltonian violates rota- 
tional symmetry. Does the violation of rotational symmetry occur also in 2+1 dimensions? 
If so, how does it manifest itself? 

In this work we initiate a systematic study of light-front QCD in 2+1 dimensions to 
investigate the various issues discussed above. The plan of the rest of this paper is as 
follows: In Sec. II we present the canonical Hamiltonian of 2+1 dimensional QCD. The 
Bloch effective Hamiltonian in the qq sector in the lowest non-trivial order is derived in Sec. 
Ill and the bound state equation is derived. The divergence structure is discussed in detail 
in Sec. IV. In Sec. V we investigate numerically, cancellation of light-front linear infrared 
divergences and the consequences of the vanishing energy denominator problem which leads 
to uncanceled infrared divergences in the bound state equation. A model which is relativistic, 
free from infrared divergences, and exhibits logarithmic confinement is presented in Sec. VI. 
In Sec. VII we present the numerical investigation of this model in the weak coupling 
limit. In this section, we also discuss the violation of rotational symmetry in this model 
at strong coupling. Finally Sec. VIII contains discussion and conclusions. Since Bloch 
effective Hamiltonian is unfamiliar to most of the readers, we present a detailed derivation 
in Appendix A. Details of the numerical procedures used in this work are given in Appendix 
B. 

II. CANONICAL HAMILTONIAN 

In this section we present the canonical light front Hamiltonian of 2+1 dimensional QCD. 
The Lagrangian density is given by 
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We have the equations of motion, 

'ij^ + gj^-m^O, (2.3) 

d^ va + gf abc A^ b Fr + g$YT a ^= 0. (2.4) 

Because we are in 2+1 dimensions, we immediately face an ambiguity since there are no 
7 matrices in 2+1 dimensions. In the literature both two component [TB[] and four com- 



ponent representation |14| have been in use. For simplicity, we pick the two component 



representation. Explicitly, 
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where £ and rj are one component fields. We choose the light front gauge A +a = 0. From 
the equation of motion, we get the equation of constraint 



id + tp~ = la 1 (id 1 + gA 1 ) + 7°m|^ 



(2.9) 



Thus the fermion constrained field 

1 
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- (id 1 + gA 1 ) + im]£. 



(2.10) 



From the equation of motion, in the gauge A +a = 0, we have the equation of constraint 
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-d x d + A la - gf abc A lb d + A lc - 2g£}T a i. 
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Using the equations of constraint, we eliminate ip and A in favor of dynamical field ip~ 
and A 1 and arrive at the canonical Hamiltonian given by 
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H = H + H int = f dx'dx^Ho + H 
The free Hamiltonian density is given by 

n = e t -(^) 2 + ^ 2 e + ^A^d'A 1 *. 

io + 2 

The interaction Hamiltonian density is given by 

Hint — T~ti + 7~(-2 
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and 



+ 2g 2 r bc ^T a O^(A lb d + A lc ) 

+ -g 2 f abc f ade ^-{A lb d + A lc )^-(A ld d + A le ). 



The one component fermion field is given by 

dk + dk l 
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b(k)e- ik - x + cf(k)t 



2(2ir)Wk+ 
The Fock operators obey the anti commutation relation 

{b(k), b\q)} = 2(2n) 2 k + 5 2 (k - q), {d(k), d\q)} = 2(2ir) 2 k + 5 2 (k - q), 



(2.16) 



(2.17) 
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other anti commutators being zero. Note that in two component representation, light front 
fermions do not carry helicity in 2+1 dimensions. 

In free field theory, the equation of motion of the dynamical field A 1 is the same as that 
of a free massless scalar field [ff5| and hence we can write 



A\x + = 0,x",x 1 ) 



dk+dk 1 



a(k)e- lk - x + a\k)e 



-ik-x 
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The Fock operators obey the commutation relation 

[a(k),a\q)] = 2(2ir) 2 k + 5 2 (k-q) 



(2.19) 



(2.20) 



other commutators being zero. 

We substitute the Fock expansions, Eqs. ( p. 17] ) and ( p,19| ) into the Hamiltonian and 
treat all the terms to be normal ordered. Thus we arrive at the canonical Hamiltonian in 
the Fock basis. 



III. BLOCH EFFECTIVE HAMILTONIAN IN THE MESON SECTOR AND THE 

BOUND STATE EQUATION 

In this section we evaluate the Block effective Hamiltonian to the lowest non-trivial order 
for a meson state and derive the effective bound state equation. We define the P space to be 
qq sector of the Fock space and Q space to be the rest of the space. In the lowest non-trivial 
order, the Bloch effective Hamiltonian is given by (see Appendix A for details) 



1 



(i | H eff | j) = (i | (# + H int ) | j) + - X> I v | k)(k | v | j) 
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1 
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(3.1) 



The states | i) and | j) are, explicitly, 
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b)=b\p 3 ,P)d\p 4 ,[3)\0), 



(3.2) 



where p±, p 2 denote momenta and a, (3 denote color which is summed over. Explicitly, 
Pi = (ptiPi) e ^ c -5 where pf is the plus component and p\ is the transverse component. For 
simplicity of notation, we will denote the transverse component of momenta without the 
superscript 1. 

The free part of the Hamiltonian leads to the matrix element 
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From the four fermion interaction, we get the contribution 
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(3.4) 



Next we evaluate the contribution from the second order term. The intermediate state 
| k) is a quark, anti-quark, gluon state. This intermediate state gives rise to both self energy 
and gluon exchange contributions. 

The self energy contributions are 

g 2 C f cW pt2(2n) 2 5 2 (p 1 -p 3 ) pj2(27r) 2 5 2 (p 2 -p 4 ) 
dkfdki 
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with 
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The gluon exchange contributions are 
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^(pt-pi) \ 1 + 1 1 , 37) 
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I pj (P^-Pi") Pj Pj (P^-pf) Pi ) 

After the construction of H e ff in the two particle space, we proceed as follows. Consider 
the bound state equation 

M 2 + P 2 
Heff I *> = p+ I *> (3.8) 

where P + , P, and M are the longitudinal momentum, the transverse momentum and the 
invariant mass of the state respectively. The state | \I/) is given by 

l*>=E f-^= f^±=MP;P,PA)bKp,mPA,P)\o) 

p J y/2(27r)*pt J y/2(2ir)*pt 



l 2(2n) 2 P+5 2 (P-pz-p A ) (3.9) 

which we symbolically represent as 

I*> = E<Mj>. (3-10) 

3 

Taking projection with the state (i \— (0 ] d(p2, a)b(pi, a), we get the effective bound state 
equation, 

M 2 + P 2 

— 02i = H 0i (p2i + Y^{l I Hi ef f I j) <j) 2 j- (3.11) 

3 

Introduce the internal momentum variables (x,k) and (y,q) via pi = xP + , p\ = xP + k, 
P2 = (l-x)P + , p 2 = (l-x)P-k, pt = yP + , p 3 = yP+q, pi = (l-y)P + , pa = (l-y)P-q 
and the amplitude 4>2(P',Pi,P2) = Trs+faix , k) . 

The fermion momentum fractions x and y range from to 1. To handle end point singu- 
larities, we introduce the cutoff rj < x,y < 1. This does not prevent the gluon longitudinal 
momentum fraction x — y from becoming zero and we introduce the regulator 5 such that 
I x — y |> 5. To regulate ultraviolet divergences, we introduce the cutoff A on the relative 
transverse momenta k and q. We remind the reader that in the superrenormalizable field 
theory under study, only ultraviolet divergence is in the fermion self energy contribution 
which we remove by a counterterm before discretization. 

The bound state equation is 
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The boson exchange contribution 
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(3.12) 
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IV. DIVERGENCE STRUCTURE 

In this subsection we carry out a detailed analysis of the divergence structure of the 
effective bound state equation. We encounter both infrared and ultraviolet divergences. 

A. Ultraviolet Divergences 

First consider ultraviolet divergences. In the super renormalizable field theory under 
consideration, with the terms appearing in the canonical Hamiltonian as normal ordered, 
ultraviolet divergence is encountered only in the self energy contributions. To isolate the 
ultraviolet divergence, we rewrite the self energy integrals as 
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2(2vr) 2C,/ 
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(y — x) 2 (l — x)(l — y) [k(l — y) — q(l — x)] 2 + m 2 (x — y)' 



(4.1] 



The first term inside the square brackets in the above equation is ultraviolet divergent, which 
we cancel by adding an ultraviolet counterterm given by 
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(4.2) 



After the addition of this counterterm, the bound state equation is ultraviolet finite. 



B. Infrared Divergences 

The infrared divergences that appear in the bound state equation are of two types: (1) 
light front infrared divergences that arise from the gluon longitudinal momentum fraction 
x g = 0, (2) true infrared divergences that arise from gluon transverse momentum k g = 
and gluon longitudinal momentum fraction x„ = 0. 



1. Cancellation of Light-front Infrared Divergences in the Effective Bound State Equation 

First consider light front infrared divergences. The effective bound state equation Eq. 
( |3.12j ) explicitly has a linear light front infrared divergent term j— -, 2 coming from instan- 
taneous gluon exchange. The most divergent part of the numerator of the transverse gluon 
exchange term in this equation is — 4 L~ g L - After combining the terms, the linear infrared 
divergent term is completely canceled and the resultant effective bound state equation takes 
the form 



M 2 - 



m 2 + k 2 



x(l — x) 
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g- 



x 



2(2tt) 2 
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The self energy contribution, made ultraviolet finite by the addition of the counterterm is 

■'-' f x , f+ A , 4m 2 
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2(2tt) 2 j J x a J-a ^[k(l-y)-q(l-x)] 2 + m 2 (x-y) 2 ' 
The energy denominator factors are 



(4.4) 
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J_ = xy J_ (1 -x){l-y) 

E\ [ky — qx} 2 + m 2 (x — y) 2 ' E% [k(l — y) — q(l — x)} 2 + m 2 (x — y) 2 



(4.5) 



The vertex terms are 

Vi=6(x-y)U(k,x,q 1 y) + 0(y - x)U(q,y,k,x), (4.6) 

V 2 = 6(x - y)U(k, 1 - x, q, 1 - y) + 9(y - z)tf (g, 1 - y, A;, 1 - x), (4.7) 

with 
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In addition to the -\ singularity which is canceled, transverse gluon exchange contributions 
also contain — singularity which is removed by the principal value prescription. Cancellation 

Xg 

of this singularity is an appealing feature of the Bloch effective Hamiltonian in contrast to 
the Tamm-Dancoff effective Hamiltonian where the singularity cancellation does not occur 
because of the presence of invariant mass in the energy denominator [ |i~6|1 . 

2. "True" infrared divergences 

Next we consider true infrared divergences. Consider the self energy integrals. The 
energy denominators in these expressions vanish when k = q and x = y which correspond 
to vanishing gluon momentum. By carrying out the integrals explicitly, in the limit A — > oo 
we get, 



SE1 = -^-C f 

2tt } 



mg 2 ^ r 1 x 1 1 — x 

x S 1 — x 5 



(4.9) 



Thus the singular part of self energy is 



mg 2 



SEl smgular = - -*- C f — r In 5. (4.10) 

2tc x(1 — x) 

The infrared divergent contribution from self energy gives a positive contribution to the 
fermion mass. It is important to note that the vanishing of energy denominator is possible 
also in 3+1 dimensions, but in that case we do not encounter any divergence. It is the 
peculiarity of 2+1 dimensions that the vanishing energy denominators cause a severe infrared 
divergence problem. 

The same vanishing energy denominators occur also in the one gluon exchange con- 
tributions. Let us now consider various terms in the numerator separately. The terms 
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proportional to 4m 2 arose from the denominator of the transverse gluon exchange. A 
straightforward calculation shows that this term leads to both finite and infrared divergent 
contributions. The infrared divergent contribution is given by 

— Cf 7 l s In 6 (4.11) 

which exactly cancels the infrared divergent contribution from self energy. The finite part, 
in the nonrelativistic limit, can be shown to give rise to the logarithmically confining po- 
tential. Next we have to consider the remaining terms in the numerator. Rest of the terms 
proportional to m 2 are multiplied by (x — y) 2 so that they do not lead to an infrared di- 
vergence problem. The numerator of the imaginary part vanishes at k — q, and x = y and 
hence is also infrared finite. It is easy to verify that the rest of the (transverse momentum 
dependent) terms in the numerator does not vanish when the denominator vanishes and 
hence the resulting bound state equation is inflicted with infrared divergences arising from 
the vanishing energy denominator. This problem was first noted in the context of QED in 
2+1 dimensions by Tarn, Hamer, and Yung |Tj| but was not investigated by these authors. 



We remind the reader that this is a peculiarity of 2+1 dimensions which provides us a unique 
opportunity to explore the consequences of the vanishing energy denominator problem. 

V. NUMERICAL STUDY OF THE BOUND STATE EQUATION 

We convert the integral equation into a matrix equation with the use of Gaussian Quadra- 
ture. (For details of the numerical procedure see Appendix B.) Cf is set to 1 for all the 
numerical calculations presented. As mentioned before, an important feature of gauge the- 
ories on the light-front is the presence of linear infrared divergences. They appear in the 
canonical Hamiltonian in instantaneous four fermion interaction term. When the qqg states 
are integrated out completely in perturbation theory, they also appear in the effective four 
fermion interaction and cancel against each other. Non- cancellation of this divergence is a 
major feature of similarity renormalization approach. We first address the issue of how lin- 
ear divergences manifest in the non-uniform grid of the Gaussian Quadrature and how well 
it can handle linear light front infrared divergence. We have studied numerically discretized 
versions of Eq. ( |3.12j ) where the divergences are present separately in the discretized version 
together with the counterterm given in Eq. (|4.2|) . For g = .2, we have calculated the eigen- 
values with and without the instantaneous interaction. The results presented in Fig. 1(a) 
for the lowest eigenvalue shows that the Gaussian Quadrature can handle the cancellation 
very efficiently. 

After the cancellation of linear light-front infrared divergence, a logarithmic infrared 
divergence which arises from the vanishing energy denominator survives in the bound state 
equation. Here we have to distinguish two types of terms. First type, where the coefficient 
of the logarithmic infrared divergence is independent of the fermion transverse momentum 
and the second type where the coefficient is dependent. Self energy and Coulomb interaction 
are of the first type. In the weak coupling limit, since the wavefunction is dominated by 
very low transverse momentum, we anticipate that contributions of the second type will be 
dynamically suppressed even though both are multiplied by the same coupling constant. 
This is especially true of any discrete grid which automatically imposes a lower limit on the 
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smallest longitudinal momentum fraction allowed. Thus at weak coupling, even if there are 
uncanceled infrared divergences (divergences of the second type), they may not be significant 
numerically whereas divergences of the second type are significant. By switching the self 
energy contribution off and on, we have studied this interplay. The lowest eigenvalue with 
and without self energy contribution is plotted in Fig. 1(b). This shows the cancellation 
of the dominant logarithmic infrared divergence. Since there are still uncanceled infrared 
divergences in the bound state equation (with coefficient proportional to fermion transverse 
momenta) this figure further illustrate the fact that such divergences are not numerically 
significant at weak coupling. 

As the strength of the interaction grows, wavefunction develops medium to large trans- 
verse momentum components and the infrared catastrophe triggered by the vanishing energy 
denominator becomes manifest numerically. This is illustrated in Table I where we present 
the variation with S of the first five eigenvalues for two different choices of the coupling g. 
The table clearly shows that on a discrete grid, the uncanceled infrared divergences due to 
the vanishing energy denominator problem are not numerically significant at weak coupling 
but their effect is readily felt at a stronger coupling. 

VI. REDUCED MODEL 

In this section we consider a model Hamiltonian free from infrared divergences con- 
structed by dropping the transverse momentum dependent terms from the numerator of the 
effective Hamiltonian. For convenience, we further drop the terms proportional to (x — y) 2 
and the imaginary part. This defines our reduced model which is also ultraviolet finite. The 
equation governing the model is given by 



M 2 - "f + k \ 
x(l — x) 



ij 2 (x,k)= SE1 ip 2 (x,k) + BE. (6.1) 



The self energy contribution SE1 is the same as given in Eq. ( |4.4| ). The boson exchange 
contribution BE is given by 

BE= -\2kr Cf L dy C d \ky - qx )*7m*{x - y) 2 MyA) 

1 n 2 r 1 r+ A Am 2 

" A dy dq— , — — r=Mv,<l)- ( 6 - 2 ) 



2 2(2vr) 2 J Jo J-A ^[k(l-y)-q(l-x)] 2 + m 2 (x~y) 2 

Again we discretize the Eq. fl6.1|) by Gaussian Quadrature. The convergence of the 
eigenvalues as a function of the number of grid points is presented in Table II. In this table 
we also present the (in) dependence of eigenvalues on the momentum cutoff. 2+1 dimensions 
provide an opportunity to study the manifestation and violation of rotational symmetry in 
light front field theory in a simpler setting compared to 3+1 dimensions. The absence of spin 
further facilitates this study. Rotational symmetry in this case simply implies degeneracy 
with respect to the sign of the azimuthal quantum number I. Thus we expect all I ^ states 
to be two fold degenerate. 

By a suitable change of variables, one can easily show that our reduced model, in the 
nonrelativistic limit reduces to Schroedinger equation in two space dimensions with a loga- 
rithmic confining potential. In the weak coupling limit, since Cf is set to 1, we can compare 
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our results of the reduced model (where we do not make any nonrelativistic approxima- 
tion) with the spectra obtained in nonrelativistic QED 2 +\- Tarn et al. [|H| solved the radial 



Schrodinger equation in momentum space for / = states and Koures |I9[ solved the coor- 
dinate space radial Schrodinger equation for general /. Since we are solving the light front 
bound state equation, rotational symmetry is not at all manifest. However, at weak cou- 
pling we expect that the spectra exhibit rotational symmetry to a very good approximation. 
Our numerical results are compared with those of Koures in Table III for two values of the 
coupling. At g — 0.2 we find reasonable agreement with the degeneracy in the spectrum. 
Even at g — .6 the violation of rotational symmetry is very small. Splitting of levels which 
are supposed to be degenerate become more visible at very strong coupling as can be seen 
from Table IV for g = 5. 

Along with the eigenvalues, the diagonalization process also yields wavef unctions. We 
have plotted the wavefunctions corresponding to the first four eigenvalues in Fig. 2 as a 
function of x and k. All wavefunctions are normalized to be J dx J dk ip 2 (x,k) = 1. The 
lowest state is nodeless and corresponds to / = 0. The next two states correspond to I = 1 
and have one node. It is interesting to note the way the node appears in the wavefunctions 
which correspond to degenerate levels. Since the rotational symmetry cannot be manifest in 
the variables x and k, how can the wavefunctions still indicate this? From Fig. 2, it is clear 
that the way this problem is resolved is by one wavefunction having a node in k and the 
other wavefunction having a node in x. Thus even if we did not know about the underlying 
symmetry from other means, the light-front wavefunctions have a subtle way of indicating 
the symmetry. 

VII. SUMMARY, DISCUSSION AND CONCLUSIONS 

In light-front Hamiltonian approach to the bound state problem in gauge theories, the 
Bloch effective Hamiltonian has certain advantages compared to the Tamm-Dancoff or the 
Bloch-Horowitz formalism. Furthermore, the recently proposed similarity renormalization 
approach is a modification of the Bloch approach. In order to quantitatively estimate the 
impact of similarity form factors in the similarity renormalization approach, it is extremely 
useful to have a quantitative study of the bound state problem in Bloch formalism. As far 
as we know, Bloch effective Hamiltonian has never been investigated in the context of the 
bound state problem in light-front field theory. 

To avoid complexities due to ultraviolet divergences we turn to 2+1 dimensions. This 
allows us to investigate light-front infrared divergences in the bound state problem in the 
presence of transverse dynamics without the additional complication arising from the mixing 
of ultraviolet and light-front infrared divergences. Further, 2+1 dimensions allow us to 
quantitatively study the manifestation and possible violation of rotational symmetry in 
light-front theory in a simpler setting. The emergence of a logarithmic confining interaction 
in the limit of heavy fermion masses is an added impetus to study gauge theories in 2+1 
dimensions. 

Only very recently, study of various issues that arise in the numerical computations in the 
similarity approach has begun. Since similarity renormalization approach is a modification 
of Bloch effective Hamiltonian approach, a detailed numerical study of the latter can serve 
as benchmark against which one can evaluate the merits of the similarity approach. It is also 
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important to quantitatively evaluate the strengths and weaknesses of numerical procedures 
in handling singular interactions in the context of light-front dynamics on the computer. 

In this work we have focused on the Gaussian Quadrature (GQ) which is one straight- 
forward procedure to solve the integral equation by converting it into a matrix equation. 
We have demonstrated the efficiency of the GQ method in handling linear and logarithmic 
light-front infrared divergences. 

A major advantage of the similarity approach is that it avoids the vanishing energy 
denominator problem that is present in the Bloch formalism. In 2+1 dimensions the van- 
ishing energy denominator leads to severe infrared divergences and hence we are presented 
with a unique opportunity to study its consequences. We encounter two types of infrared 
divergences: (1) with a coefficient proportional to fermion mass and (2) with a coefficient 
proportional to fermion transverse momentum. The former type gets canceled in the bound 
state equation between fermion dressing by gluon and gluon exchange between fermions. 
The latter type is uncanceled which however can be dynamically suppressed at very weak 
coupling on a finite grid. We have demonstrated that on a discrete grid provided by GQ, the 
uncanceled divergences are numerically insignificant at weak coupling whereas the catastro- 
phe due to their presence is readily felt at stronger coupling. 

We proceed to study a reduced model that is free from infrared divergences and which 
reduces to the Schroedinger equation with a logarithmic potential in the nonrelativistic 
limit. This model provides us an opportunity to study the simplest manifestation and 
possible violation of rotational symmetry in the context of light-front field theory. Even 
though the Hamiltonian does not exhibit rotational symmetry we have shown that at weak 
coupling spectra exhibit rotational symmetry to a very good approximation. We have also 
shown that even though the rotational symmetry is not manifest in the variables x and k, 
light-front wavefunctions have a subtler way of indicating the underlying symmetry. 

Our study indicates that in the context of Fock space based effective Hamiltonian meth- 
ods to tackle gauge theories in 2+1 dimensions, approaches like similarity renormalization 
method are mandatory due to uncanceled infrared divergences caused by the vanishing en- 
ergy denominator problem. It is important to recall that Bloch effective Hamiltonian is 
generated by completely integrating out the intermediate gluons irrespective of whether 
they are low energy or high energy. Is this justified in a confining theory? 

Now that we have obtained quantitative measures of the vanishing energy denominator 
problem and the nature of the spectra at weak coupling of the Bloch effective Hamiltonian, 
the next step is to study QCD 2 +i in the similarity renormalization approach which avoids 
the vanishing energy denominator problem. An important issue here is the nature of new 
effective interactions generated by the similarity approach. It has been shown that in 3+1 
dimensions, similarity approach generates logarithmic confining interactions which however 
breaks rotational symmetry. It is interesting to investigate the corresponding situation in 
2+1 dimensions. 

APPENDIX A: BLOCH PERTURBATION THEORY FOR EFFECTIVE 

HAMILTONIAN 

Bloch perturbation theory was introduced in Ref. 0. Here we follow the treatment in 
Ref. |T7] where the reader can find many examples of perturbative calculations. 
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Consider a Hamiltonian H defined at a cutoff A. Let us try to lower the cutoff to A. 
In general, the cutoff could be in energy and/or particle number. Let us denote by Q 
the operator that projects on to all of the states removed when the cutoff is lowered. Let 
P = I -Q. We have 



Q 2 = Q, P 2 = P, PQ = QP = 0. 



(Al) 



Our purpose is to find an effective Hamiltonian H e ff that produces the same eigenvalues in 
the sub space P as the original Hamiltonian H. 
Introduce an operator R that satisfies 



Q | V) = RP I V>) 



(A2) 



for all eigenstates of the Hamiltonian that have support in the subspace P. R gives the part 
of | ip) outside the space projected by P in terms of the part of | ip) inside the space. Require 
that R gives zero acting on states outside the subspace. This means R = RP, R = QR, 
R 2 = 0. From R = QR, we have, PR = 0. Note also that R ] ^ R. 
Start from the set of equations 



From Eq. 



From Eq. (K3) 



PHP | ip) + PHQ | if>)= EP | ip), 
QHP\iP)+QHQ\iP)=EQ\iP). 



RPHP | ip) + RPHQRP | ip) = ERP \ ip). 



QHP | ip) + QHQRP | V) = ERP \ ip). 



Subtracting, 



RHpp — HnnR + RHpnR — H, 



PQ 1 



l QP 



0. 



(A3) 
(A4) 



(A5) 



(A6) 



(A7) 



We have introduced the notations, PHP = Hpp and so on. Put H = h + v with [h, Q] = 0. 
Then 



Rhpp - h QQ R - v Q p + Rvpp - vqqR + Rv PQ R = 



(A8) 



which shows that R starts first order in v. 
We start from the eigenvalue equation, 



H(P + Q) \iP)=E{P + Q) \ip). 



(A9) 



i.e., 



H(P + R)P | ip) = E(P + R)P | ip). 
Multiplying from the left by (P + R<) we have, 



(A10) 
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(P + R^)H{P + R)P | V) = P(P + R r ){P + R)P | # 



(All) 



Using PR = 0, R)P = 0, (P + P f )(P + P) = P + P f P. Thus we can rewrite the eigenvalue 
equation as 



-1 + PtP 
i.e., 

where 

and 



: (P + R^H(P + R) 



1 + PtP 



'[1 + tfRfip | V) = P[l + R ] R\^P | # (A12) 



H, 



eff 



E 



H, 



eff 



1 



1 + PtP 



[l+P f P]^P|^) 



: (P + tf)H{P + R) 



1 il 



1 + PtP 



(A13) 



(A14) 



(A15) 



Our next task is to generate a perturbative expansion. Denote free eigenstates in P by 
a), | b), etc. Denote free eigenstates in Q by | i), \ j), etc. Then 



h PP | a)= e a | a), 



(A16) 



Let us compute P to lowest orders in the perturbation theory. Let us write P = P1 + P2 + . . 
where the subscript denotes orders in v. A straightforward calculation leads to 



(i I Pi I o)- 



(* I V QP I °) 



(i|P 2 |a)=-^- 

6 {e a -ei){e b -e 



v\a)(i\v\b) (i\y\j)(j \y\g) 

i (e a - ei)(e a - €j) 



(A17) 
(A18) 



Our next task is to develop a perturbation theory expansion for the effective Hamiltonian 
to a given order. 

We start from the expression for the effective Hamiltonian. Remember that Pi ~ 0(v), 
P 2 ~ 0{v 2 ). 

To order v, H e ff — PHP and hence 



(a I H eff I b) = (a \ (h + v) \ b). 
To second order in v, we have 



(A19) 



H eff = [1 - -R^R] [PHP + PHR + R^HP + R ] HR] [1 - -R ] R\. (A20) 



From R^HR we get, 
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{a\n?HR\Q = 'E« V VU \? ] "'? ■ < A21 » 

■ (e„-ei)(e 6 - ei ) 

From PHR and PJHP terms we get 

J2(a \H\i)(i\R 1 \b) + J2( a \ R l\i)(i\ H \ h ) ( A22 ) 

i i 

= y r (a | v \i)(i \ v \ b) + (a \v \i)(i \ v \ b) 



From the normalization factors we get 

- Wtf P - IpEPHfR = -i( £ „ + *) £ ^'"^^"'y (A24) 

II I i (e a — ei){e b — €i) 

Adding everything, to second order, we have, 

(« I H eff | 6) = ij> | v | »>(< | v | 6>f-^— + -^— 1. (A25) 



2 



^a tj 65 fcj 



If a = 6, this expression reduces to the familiar second order energy shift. 

Why Bloch formalism is preferred over Bloch-Horowitz formalism? 

In the former, eigenstates of the effective Hamiltonian are ortho normalized projections 
of the original eigenstates. In the latter, they are not. 

Consider two ortho normalized eigenstates of the original Hamiltonian | ipi) and | 1P2) 
with (0i I ip 2 ) = 0. However, P \ tpi) and P \ 1P2) need not be orthogonal, i.e., {ipi | PP \ 
^2) — (V'l I P I 4 ) 2) 7^ 0. Consider 

(V'i I ^2)= (in\P\ ^2) + (V'i I Q 2 1 ^2) 

= (^1 I P I V2) + (V>i I P f P f PP I Vi)- (A26) 

Construct | &) = [1 + fltfljip | ^ | ^ 2 ) = [1 + ptpjlp | ^ 2 ). Then 

$1 I ^2) = (Vi I P I V> 2 ) + (V>i I PR ] RP I ^2) = (V>i I ^2). (A27) 



APPENDIX B: DETAILS OF NUMERICAL PROCEDURE 

Parametrization: The light-front variables are parametrized in the following ways in our 
numerical calculations. The full /c-interval is divided into n\ quadrature points, k is defined 
by two different ways. One definition is 

uAm (Bl) 



;i-M 2 )A + m' 



where A is the ultraviolet cutoff and w's are the quadrature points lying between —1 and 
+ 1, so that k goes from —A to +A. The other definition is 



1 7/7T 

k = -tan(^), (B2) 

here k is a parameter that can be tuned to adjust the ultraviolet cutoff. The second definition 
( |B2| ) of k is very suitable for weak coupling calculations where we need maximum points to 
be concentrated near k = and get better convergence than the first definition (|Bf| ). 

The longitudinal momentum fraction x ranges from to 1. We divide all x- integrations 
in our calculations into two parts, x ranging from to 0.5 and x ranging from 0.5 to 1 and 
discretize each ^-interval into n2 quadrature points with the parametrization 

x = 1 + " + 2 J (1 ' t;) , n < x < 0.5, (B3) 

o;= 3 + t; '^ (1 + " ) , 0.5<x<l-^ (B4) 

where v's are the Gauss-quadrature points lying between —1 and +1 and r/(— > 0) is intro- 
duced to handle end-point singularities in x as mentioned before. 

To handle the infrared diverging terms we put the cutoff \x — y\ > 8 and at the end we 
take the limit 5 — > 0. Numerically, it means that the result should converge as one decreases 
5 if there is no net infrared divergence in the theory. 

Diagonalization: After discretisation, solving the integral equation becomes a matrix 
diagonalization problem. The diagonalization has been performed by using the packed 
storage LAPACK p0[ routines DSPEVX for the reduced model (real symmetric matrix) 
and ZHPEVX for the full Hamiltonian (Hermitian matrix). 
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FIGURES 

FIG. 1. Cancellation of infrared divergence. Full line denotes the full Hamiltonian. (a) shows 
the cancellation of light-front infrared divergence by switching on and off the instantaneous 
interaction. Filled circles - without instantaneous interaction, (b) shows the cancellation of 
logarithmic infrared divergence by switching on and off the self energy term. Filled circles 
- without self energy. The parameters are g = .2, rj = .00001, m — 1, K — 20, ri\ = 40, 
n 2 = 50. 



FIG. 2. The wavefunctions corresponding to the lowest four eigenvalues of the reduced 
model as a function of x and k. The parameters are g — .2, r] — .00001, m — 1, k — 10, 
rii = 46, n 2 = 74. (a) Lowest state, (b) first excited state, (c) second excited state, (d) third 
excited state. The first and second excited states should be degenerate in the absence of 
violation of rotational symmetry. 
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9 


5 


eigenvalues (M 2 ) 


0.2 


0.00001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.0001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.005 


4.0901 


4.1066 


4.1099 


4.1100 


4.1112 


0.01 


4.0870 


4.0972 


4.0972 


4.0973 


4.0973 


0.6 


0.0001 


-187230.4 


-187225.4 


-186664.9 


-186664.8 


-31506.9 


0.001 


-187230.4 


-187225.4 


-186664.9 


-186664.8 


-31506.9 


0.005 


1.9094 


1.9415 


3.1393 


3.1399 


4.5697 


0.01 


4.5735 


4.7337 


4.7667 


4.7832 


4.8277 



TABLE I: Variation with 5 of the full hamiltonian. The parameters are nl=40, n2=50, 
77=0.00001, k=20.0 ink= \tan(^f) 
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Ill 


n2 


eigenvalues (lowest five) (k =10.0) 


20 


20 


4.08926 


4.10605 


4.10768 


4.11061 


4.11085 


30 


30 


4.09045 


4.10909 


4.11038 


4.11516 


4.11699 


40 


30 


4.09045 


4.10913 


4.11035 


4.11524 


4.11697 


40 


40 


4.09102 


4.11052 


4.11154 


4.11711 


4.11951 


40 


50 


4.09136 


4.11133 


4.11222 


4.11811 


4.12096 


50 


50 


4.09136 


4.11135 


4.11219 


4.11816 


4.12095 


50 


60 


4.09158 


4.11188 


4.11263 


4.12189 


4.12290 


46 


60 


4.09158 


4.11187 


4.11264 


4.11877 


4.12189 


46 


66 


4.09168 


4.11212 


4.11284 


4.11905 


4.12231 


46 


74 


4.09179 


4.11237 


4.11305 


4.11934 


4.12276 


nl 


n2 


eigenvalues (lowest five) (A=20.0) 


46 


74 


4.09179 


4.11240 


4.11301 


4.11940 


4.12273 



TABLE II: Convergence of eigenvalue with nl and n2 (reduced model). The parameters 
are m=1.0, g=0.2, rj = 0.00001. 
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9 


eigenvalues 


0.2 


This 

work 


4.0918 
(4.1227, 4.1235) 


(4.1124, 4.1130) 
(4.1268, 4.1273) 


4.1194 
(4.1298, 4.1303) 


Koures 
(Ref. 0) 


4.0925 (1 = 0) 
4.1260 (I = 2) 


4.1144 (I = 1) 
4.1303 (I = 1) 


4.1214 (I = 0) 
4.1340 (1 = 3) 


0.6 


This 

work 


4.5856 
(4.8767, 4.8816) 


(4.7741, 4.7821) 
(4.9094, 4.9184) 


4.8390 
(4.9458, 4.9481) 


Koures 
(Ref. H) 


4.5806 {I = 0) 
4.8827 (/ = 2) 


4.7777 (/ = 1) 
4.9205 {I = 1) 


4.8409 (Z = 0) 
4.9545 {I = 3) 



TABLE III: Reduced model. The parameters are nl=46, n2=74, rj = 0.00001, m=1.0. 

k = tan(q7r/2)/K, k = 20.0. Eigenvalues within () are ±Z degenerate (broken) states. 
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111 


n2 


eigenvalues 


I 


40 


50 


18.217 


(30.702, 33.499) 


35.206 


(39.955, 41.159) 


(41.332, 43.271) 


(44.134, 45.272) 


46 


70 


18.276 


(30.774, 33.616) 


35.318 


(40.106, 41.331) 


(41.483, 43.477) 


(44.375, 45.503) 


II 


40 


50 


18.980 


(31.507, 34.219) 


35.826 


(40.406, 41.888) 


(41.921, 43.788) 


(44.345, 45.163) 


46 


70 


19.008 


(31.542, 34.319) 


35.935 


(40.626, 42.031) 


(42.088, 44.010) 


(44.647, 45.780) 



TABLE IV: First few eigenvalues in the reduced model. The parameters are g=5.0, 
m=1.0, r/=0.00001. (I) for the parametrization k = uAm/((l - u 2 )A + m), A = 40.0. (II) 
for the parametrization k = tan(u7r/2)/K, k = 10.0. 
Eigenvalues within ( ) are ±/ degenerate (broken) states. 
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FIG. 1. Cancellation of infrared divergence. Full line denotes the full Hamiltonian. (a) shows 
the cancellation of light-front infrared divergence by switching on and off the instantaneous inter- 
action. Filled circles - without instantaneous interaction, (b) shows the cancellation of logarithmic 
infrared divergence by switching on and off the self energy term. Filled circles - without self energy. 
The parameters are g = .2, 77 = .00001, m = 1, k = 20, rii = 40, ri2 = 50. 
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FIG. 2. The wavefunctions corresponding to the lowest four eigenvalues of the reduced model 
as a function of x and k. The parameters are g = .2, rj = .00001, m = 1, k = 10, n\ = 46, n^ = 74. 
(a) Lowest state, (b) first excited state, (c) second excited state, (d) third excited state. The first 
and second excited states should be degenerate in the absence of violation of rotational symmetry. 
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